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Abstract 

Decay of regular static spherically symmetric solutions in the SU{2) Yang-Mills-dilaton (YMd) 
system of equations under the independent excitation of their unstable eigenmodes has been stud- 
ied self-consistently in the nonlinear regime. The considered regular YMd solutions form a discrete 
family and can be parametrised by the number N = 1,2,3,4... of their unstable eigenmodes 
in linear approximation. We have obtained strong numerical evidences in favour of the follow- 
ing statements: i) all static YMd solutions arc distinct local threshold configurations, separating 
blowup and scattering solutions; ii) the main unstable eigenmodes are only those responsible for 
the blowup/scattering alternative; iii) excitation of higher unstable eigenmodes always leads to 
finite-time blowup; iv) the decay of the lowest N = 1 static YMd solution via excitation of its 
unique unstable mode is an exceptional case because the resulting waves propagate as a whole 
without energy dispersion revealing features peculiar to solitons. Applications of the obtained 
results to Type-I gravitational collapse of massless fields are briefly discussed. 

PACS numbers: 1.27.+d; 03.65.Pm; 02.60.Lj 



I. INTRODUCTION. 



In the last decade a considerable progress has been achieved in understanding of singular 
"blowup" solutions of nonlinear evolution equations with smooth initial data possessing finite 
energy. Some issues under consideration were inspired by pioneering work by Choptuik [1] 
on a gravitational collapse of massless field. This task has been developed in a number of 
new ways including gravitational collapse of nonlinearly interacting massless fields, such as 
Yang-Mills field and nonlinear sigma-models [2] , [3] , [1] . Other related issues were motivated 
by recent studies on blowup in nonlinear wave equations [3] - |16j . 

The nature of the threshold of singularity formation is of particular interest. This problem 
was studied analytically and, mostly, numerically and it was realized that there are at least 
two kinds of critical solutions which play a role of local thresholds, separating dispersing 
and singular solutions. 

The thresholds of first kind are self-similar solutions possessing one unstable mode. Criti- 
cal solutions of this type can be both discretely self-similar or continuously self-similar. Con- 
tinuously self-similar thresholds are well-known solutions and their explicit form was found 
in some cases even analytically as solutions of system of ordinary differential equations after 
substitution of the "self-similar" anzats. Discretely self-similar critical solutions are more 
subtle objects and they were discovered first in numerical studies of Type-II gravitational 
collapse of massless fields [T], [S], [1], [T7]. Black holes of arbitrary small mass (Type-II 
gravitational collapse) can be produced in this way and the critical discretely/continuously 
self-similar solution arises as a threshold, separating formation of naked singularity (zero- 
mass black hole) from dispersing solutions. 

Another kind of the threshold configurations are static solutions with at least one unstable 
eigenmode. The well-known example of such a threshold is lowest = 1 Bartnick-McKinnon 
(BK) solution in a coupled system of Einstein- Yang-Mills equations [IS]. This threshold 
solution provides Type-I gravitational collapse scenario with a finite mass gap in a spectra 
of produced black holes; the smallest black hole, formed in this way has a mass equal to 

= 1 BK static solution [2]. 

In this paper we focus on static threshold configurations. We consider spherically sym- 
metric system of coupled SU (2) Yang-Mills-dilaton (YMd) equations in 3 + 1 Minkowski 
space-time. This system possess infinite number of finite energy static regular solutions 
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which are parameterized by the number of zeros = 1,2,3,4... of relevant YM function 
[T^ . All of these static solutions are unstable and have N even-parity unstable eigenmode 
in a small linear perturbations approach. We have shown in one of our previous papers 
[9] that this YMd system is supercritical in PDE terminology ^20j and all solutions of the 
corresponding evolutional Cauchy problem with "large enough" initial data result in finite 
time blowup. 

It is widely believed that a static unstable solution with one unstable eigenmode could 
be a threshold configurations in supercritical PDE equations. In other words, these static 
solutions should have stable manifold of codimension-one and the instability results from 
a simple negative eigenmode. This conjecture was proven in some cases analytically [H] 
and numerically [I5]. It was shown numerically in various models that blowup/scattering 
alternative in a final solutions asymptotics results from the sign of the excited single unstable 
eigenmode [15], [21] . 

The considered YMd system represents more general case because it posessess an infinite 
family of the static solutions with = 1, 2, 3, 4, . . . oo unstable modes correspondingly. A 
stable manifold of these solutions with A^ > 1 has a codimension more than one and few 
natural questions arise as a result: what is a role of higher static solutions with A^ > 1; do 
they also relate to the distinct thresholds, separating blowup and scattering; and if so, what 
is the role of higher unstable modes in this blowup/scattering alternative? 

To answer these questions, we systematically studied the decay of various lowest static 
YMd solutions under the independent excitations of all of their proper unstable eigenmodes. 
The results of this numerical studies are presented in Section 4. Section 2 contains basic 
notations and main field equations. In Section 3, we discuss the basic properties of the static 
intermediate attractors. Conclusions, discussion, and further open questions are summarized 
in Section 5. 

II. PRELIMINARY BASIC FACTS AND EQUATIONS. 

A coupled system of Yang-Mills-dilaton (YMd) fields in 3 + 1 Minkowski space-time 
is given by the action: 

^ = il (^(^^^' - ^^^"^^^;^) d'xdt. (1) 
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Here $ is the dilaton field, is the Yang-Mills field, and k and g are dilaton and gauge 
coupling constants, respectively. 

In the spherically symmetric case the dilaton field and the SU{2) Yang-Mills potential 
can be described in terms of two independent functions 

k 

<l> = $(t,r), A» = 0, A^ = e,,fc^(/(t,r)-l). (2) 

After substitution ^ into Q and rescaling $ ^ $/A;, r ^ {k/g)r, t — > {k/g)t, S gkS 
the dependence on two parameters k and g effectively vanishes, and the integration over the 
angular variables gives finally the reduced action: 



S 



1 1 / ('f2_1^2 



dr dt. (3) 



2" 2" ' ' y ' 2r2 

The corresponding system of field equations has the form of a coupled nonlinear wave equa- 
tions in the following way: 

f(l - P) 

ftt + frr -fr^r= \J ' , (4a) 



^U-^rr-'-^ = --Af^--f^+'-^^]- (4b) 



r2 ■'^ 2r2 

Invariance of the field equations Q with respect to the replacement / — / allows to 
put Yang-Mills function f{t, 0) = 1 at the origin without loss of generality. We look for finite 
energy solutions of Q which are regular both at the origin and at the spatial infinity. These 
requirements provide boundary conditions for the functions r) and f{t,r) compatible 
with the field equations as follows: 

/(t,r = 0) = l, /,(t,r = 0) = 0, $(t = 0,r = 0) = $o, $,(t,r = 0) = 0; (5) 
lim /(t,r) = ±1, lim /,(t,r) = 0, 

r— >oo r— >oo 

lim r) = $oo, lim r) = 0. (6) 

Here $o is an arbitrary finite constant and $oo is a constant, depending on $o and on the 
total mass-energy of the configuration. Note that according to the regularity requirement 
Yang-Mills function should take its vacuum values / = ±1 at the origin and at the spatial 
infinity. 

The total mass-energy E of the configuration described by the field equations Q is given 
by the energy functional: 

+ 00 +00 

^1 _ _ 1 _ _ - / . ^ (f2 _ T^2M 

dr, (7) 



+ 00 +00 







where the energy density function e{t,r) = Tl{t^r) is introduced in a standard way as the 
tt component of the energy-momentum tensor T^. Energy E is unique on-shell conserved 
value for the system Q: dE/dt = 0. 

For the further sake it is useful to introduce radial momentum density function p{t, r) 
which is a tr component of the energy-momentum tensor; in terms of the relevant functions 
r) and f{t,r) it has form: 

p{t, r) = T[ = r2$,(t, r)<!>r{t, r) + 2e* ft{t, r)/,(t, r). (8) 

+00 

The total radial momentum = J T[{t,r)dr is not conserved for the given system, 



since the action functional ^ contains explicitly radial variable r. However, at any given 
values of t and r this function p{t, r) describes the total radial energy flux through the 
surface of sphere of radius r and the sign of function p{t, r) indicates the direction of energy 
propagation: positive sign of p(t, r) corresponds to the ingoing (towards to the origin) energy 
flux and negative - to the outgoing energy flux. 

Static {t -independent) regular solutions of system Q with the boundary conditions (|5]), 
(|6| are well-known since early 90th [T9]. Similar to the Bartnick-McKinnon (BK) solutions 
of the SU(2) Einstein- Yang-Mills equations system [TH], there is an infinite discrete set of 
the regular static solutions {fN{^)i ^N{f)} labeled by the total number of zeros of YM 
function /(r), = 1,2,3,4, ...cxd, which were found numerically. The Yang-Mills functions 
/iv(r) oscillate — 1 times, being bounded inside the strip /7v(^) £ 1]) whereas the 
dilaton functions $Ar(r) are monotonically growing. 

The scale transformation: 

/^/, $-^$ + A, r^rexp[A/2], t ^ texp[A/2], A = const (9) 

does not change the field equations either. That is why without loss of generality one can 
put $0 = <l>(t = 0, r = 0) = in Q. The mass-energy takes a scale-dependent value and it 
transformed under the scale transformations (|9| as: 

M exp[A/2]M. (10) 

If the scale is chosen according to a lirrir^oo^Nif^) ~^ 0, then the mass-energy of the cor- 
responding static YMd solutions increases with as follows: Mi = 1.60751, M2 = 1.93112, 
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TABLE I: Unstable eigenvalues {A^}^;^ for the background static YMd solutions with N = 



1,2,3,4. 



N 










1 


-9.0566 X 10-2 








2 


-7.5382 X 10-2 


-2.0742 X 10-^ 






3 


-4.9346 X 10-2 


-1.4957 X 10-^ 


-1.9622 X 10"^ 




4 


-4.3455 X 10-2 


-5.9905 X 10-5 


-1.3278 X 10-^ 


~ -10-9 



M3 = 1.98863, M4 = 1.99814, M5 = 1.99970, Mg = 1.99995, ...,Moo = 2.0, saturating their 
upper limit Moo = 2.0 as ^ cxo [19]. 

In terms of the hnear perturbation analysis every such solution with zeros of YM 
function {/Ar(r), $Ar(r)} has even-parity unstable modes [I9]. Indeed, assume spherically 
symmetric perturbations of the regular static solutions {/jv(r), ^fyf{r)} of the following form: 

f{t,r) = fN{r) + ee-^''^'-^/^v{r)e''^\ <l>(t, r) = <l>^(r) + e — ^(r) e'"*, (11) 

r 

where e is a small parameter. This gives a standard matrix (2x2) Sturm-Liouville eigen- 
value problem with respect to the appropriately normalized eigenfunctions {u{r),v{r)}. The 
negative eigenvalues ct;2 = A < correspond to the unstable eigenmodes. All the unstable 
eigenvalues A]^, j = 1, . . . , A^ that correspond to the background solutions with A^ = 1, 2, 3, 4 
are presented in the following Table 1 in accordance to our paper |22]. The corresponding 
eigenfunctions were also found and plotted in [22], [23] and we use them to determine the 
initial data for the Cauchy problem in Section 4. 

Now we proceed to another story inspired by the last decade developments in understand- 
ing of singularity formation for several classes of nonlinear wave equations. It was shown 
in our previous paper [9] that the considered system of Yang-Mills-dilaton equations ^ 
admits a hidden scale-invariant form. This observation allows to make definite predictions 
about possible late-time solution asymptotics in a nonlinear regime under the assumption 
of a well-posed evolutional Cauchy problem. 

Indeed, if we decompose the dilaton function $(t, r) in the following way: 

$(t,r) = 0(t,r) + 21nr, (12) 
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the system of equations Q rewritten in terms of functions f{t,r) and 4>{t^r) becomes in- 
variant under the global scale transformations of time and radial coordinate t t/X , 
r r/\. The mass-energy functional ([T]) expressed in terms of the scale-invariant functions 
becomes 

E = |rfr(V0'2 + 2r0' + 2 + V02 + rV(/'2 + /2 + ^^^?^ (13) 

providing the homogeneous scale law for the mass-energy as 

t r , . ,t r 



E 



A"E[/(t,r),0(t,r)], a = +l. (14) 



According to general expectations [20] the degree a of scale parameter A, occurring in the 
energy homogeneous scale law, defines the criticality class of the considered PDE equation, 
or system of equations. In particular, if a > the system is called supercritical as this 
takes place in our case and one should expect singularity formation in the corresponding 
well-posed Cauchy problem in a finite time and for all initial data, exceeding some threshold 
value. One should also expect that a part of the solution shrinking to the origin evolves in 
certain universal way prior to the singularity formation. It is believed that the mentioned 
universality is inspired by the existence of an attractor which these shrinking solutions attain 
for a while. This attractors, as a rule, are described by self-similar solutions. Exploring 
scale properties of the equations allows to search for an explicit form of possible self-similar 
attractors. 

In our case the revealed invariance of the solutions under the scale dilations 

allows to search for the solutions which effectively depend only on a particular combination 
of the independent variables - t/r, i.e. for self-similar solutions as follows: 

/(t,r) = /(x), 0(t,r)=0(x), a; = — . (15) 



The function (f){t,r) defined in (12) and the positive constant T is transformed similarly to 
r, t under the dilations. The constant T could have the meaning of a blowup time - absolute 
value of the time in the evolution Cauchy problem when the singularity at the origin start 
to develop and the solution becomes nondifferentiable. 

In terms of the functions /(x) and 0(x) the system of PDEs ^ is transformed into a 
system of ordinary differential equations (ODEs). This system of ODEs with the appropriate 



boundary conditions was considered in [9] and the infinite discrete family of the desired 
self-similar solutions was found. These self-similar solutions can be parametrized by = 
0, 1, 2, 3,4, ... + oo - the number of zeros of the Yang-Mills function f{x). According to the 
performed linear perturbation analysis only the lowest self-similar solution with = is 
stable. 

And just this the lowest A^ = self-similar solution can pretend to be a global stable 
attractor in the evolutional Cauchy problem as was shown in [9|. The Cauchy problem 
for the system (|4]), ([5|, ([6]) has been studied numerically for a wide range of smooth finite 
energy initial data. It turns out that if the initial data exceeds some threshold, the result- 
ing solutions in a compact region that shrinks to the origin attain the A^ = self-similar 
solutions. These solutions, being self-similar, further evolve in the universal way in a finite 
time until the second derivative of the YM function at the origin starts growing infinitely 
what indicates on a singularity formation. 

It should be noted that such a singularity formation (blowup) at the origin is not ac- 



companied by energy concentration according to the energy scale law (14). This scale law 
describes the energy of a self-similar part of the solution as it evolves, shrinking to the origin, 
whereas the total energy of the whole configuration ([T]) is conserved. Indeed, the energy of 
self-similar part of the solution at time t <T {T -blowup time) inside the past light cone of 
the point t = T, r = is equal to M = E 4(T — t) and vanishes as t ^ T [9]. 

The system dynamics determined by the near-threshold initial data in the mentioned 
Cauchy problem is a task of particular interest. In fact, solutions with various initial data 
approaching their corresponding threshold values could attain for a while (prior to attain 
A^ = self-similar attractor) some universal threshold configuration. In our case the A^ = 1 
static YMd solution is one of such threshold configurations that plays a role of an interme- 
diate attractor. We discuss this in more details in the next section. 



III. NOTES ABOUT STATIC INTERMEDIATE ATTRACTORS. 

This section, as well as the previous one, contains rather preparatory material to the next 
section. We would like to show that the A^ = 1 static YMd solution of the system (|4]), ([s]), ([6]) 
is a local intermediate attractor and it is also a threshold configuration separating blowup 
and scattering solutions. Following the conventional methodology we assume an isolated 
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compact distribution of matter (YMd) fields as an initial data in the Cauchy problem (|4]), 
(|5|, Further nonlinear evolution of the corresponding initially ingoing wave is a main 
subject of present section. 

In this section we consider a compact distribution of YM field, determined by the function 
f{t = 0,r) and by their partial derivatives at the initial slice t = 0. Then the initial profile 



of the dilaton function ^{t = 0,r) is obtained by integration of the dilaton equation (4b) 
from the origin up to infinity assuming $(t,r = 0) = 0, $(t = 0,r) j = 0, $(t = 0,r) = 
(see [5] for more details ). The resulting value of the dilaton function at spatial infinity 
r oo) = <|)?^°/*'^ is a key parameter which defines the position of the virtual static 
YMd solutions in this scale. Indeed, in order to study possible attaining static solutions one 
should plot all their virtual images at the same scale, i.e. with $Af('" oo) = ^p^"-^*'"^ by 
using an appropriate shift of the dilaton function and the scale transformations ([9]). This 
was done for all simulations discussed below. 

Without loss of generahty we assume two types of the initial YM profile described by the 
one-parametric kink-type and Gauss-type distributions as follows: 

fit = 0, r) = 1^^, f{t = 0,r) = l-A exp [-a{r - R)'] , (16) 

where a is fixed and free parameters a, A control the strength of the given initial data. 

If the initial data are less than the threshold values in terms of acr, Acr, than the initially 
ingoing waves get smoothly bounced near the origin and are radiated away similarly to 
the solutions of linear wave equations. Positions of various YM and dilaton profiles at 
the moment of bouncing are plotted in the Figure 1 together with virtual images of the 
lowest = 1, 2 static YMd solution at the appropriate scale. One can see that, as we are 
approaching the critical values of the parameters towards more strength initial data, the 
position of bounce is approaching from outside the virtual = 1 static YMd solution. 

Dynamics of two various near-threshold solutions is presented on the Figure 2. Behavior 
of the asymptotically scattered solution in a vicinity of the virtual = 1 static solution is 
presented in the left column of the Figure 2. One can see that initially purely ingoing wave 
attains the virtual static = 1 YMd solution for a while and after smooth bounce from 
it radiates away being transformed to the outgoing wave. If the relevant parameter in the 
initial profile slightly exceeds the threshold value one has a blowup as a final asymptotics. 
The behavior of the corresponding solution is presented in the right column of the Figure 2. 
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dr, (17) 







The initially purely ingoing wave attains the virtual static = 1 YMd solution for a while 
too, and then starts to collapse, attaining later on the universal = self-similar solution 
which is a blowup precursor. 

The time-dependent part of the energy functional for the considered solutions 

'1 

has a local minima at some instant of time for all solutions which attain the intermediate 
attractor. Gradual change of the initial parameters a and A in a small vicinity of their 
threshold values allows to minimize this integral over all range of the initial data parame- 
ters. The corresponding positions of the solution profiles (plotted at the moment when the 
integral has a minimum) converge to the limiting (static) threshold configuration separating 
asymptotically scattered and blowup solutions. It is not surprising that the = 1 static 
YMd solution is occurred to be a such limiting threshold configuration. This step-by-step 
shooting numerical strategy allows to approach the threshold value of the relevant parameter 
in the initial data with precision about 10~^^. 

The mass-energy of the static intermediate attractor is an absolute lower bound for the 
mass-energy of all solutions of the ingoing wave type which could attain this attractor. This 



(unreachable) limit corresponds to the vanishing integral (17). Hence, in order to have a 
chance to attain a static intermediate attractor, the mass-energy of the initial field configu- 
ration must exceed the mass-energy of the corresponding virtual static solution (calculated 
at the appropriated scale determined by the <|)P^°/*''^) which pretends to be an attractor. In 
this case a considerable part of the energy surplus should radiate away prior the solution 
attains the static attractor. The observed formation of an outgoing energy flux could serve 
as a useful numerical tool to resolve a possible attaindance to the intermediate attractor. 
The bottom plots of the Figure 2 illustrate this statement: the negatively deflned parts of 
the radial momentum density function p{t,r) testify to occurrence of the outgoing energy 
fluxes which are formed prior the solutions attain the static = 1 intermediate attractor. 

To summarize: at least for those solutions which temporary attain the virtual A^ = 1 
static YMd solution the last one plays a role of a local threshold conflguration separating 
scattering and blowup solutions. This is an expected behavior for the supercritical systems 
admitting static solutions with one unstable eigenmode. Note that this is the most clear 
case, since due to an evident transversality, a stable manifold of the A^ = 1 static solution 
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has codimension 1. Under the circumstances it was possible to observe convergence to the 
= 1 static threshold configuration starting with the arbitrary chosen simple kink-type or 
Gauss-type initial data at spatial infinity. 

Situation becomes more complicated if we consider higher static YMd solutions with 
> 1 which have more than one unstable eigenmodes. Stable manifold of these solutions 
with the > 1 has codimension > 1 as well and in this case, for a number of reasons, it 
is difficult to create initial data at spatial infinity which could control (by the fine tuning 
of the governing parameters) the possible attaindance to the higher > 1 static solutions 
pretending to be distinct threshold configurations. 

Instead, it is reasonable to start directly with this unstable static solutions and study 
their decay via unstable eigenmodes excitation. The role of various negative eigenmodes 
studied from the viewpoint of blowup/scattering alternative is a main task of this paper and 
is considered in the next section. 



IV. DECAY OF STATIC YANG-MILLS-DILATON SOLUTIONS AND THEIR 
ASYMPTOTICS. 

We consider the coupled system of the YMd nonlinear wave equations Q with the 
boundary conditions (|5|, The initial data at t = are determined by the static 

solution {/jv(r), $Ar(r)} (A^ is fixed) perturbed by their proper unstable eigenfunctions 
{f]v(r), M5y(r)}, {i is also fixed, 1 < i < A^) with the corresponding eigenvalue A^y = < 



according to (11) and Table 1 as follows: 



fit = 0,r) = /^(r) + ee-*-M/2^ivW, Ht = 0,r) = $^(r) + e ^ nV(r) , (18) 

r 

where e is small parameter. Time derivatives at the initial time slice t = are obtained 



from the (|TT|) as follows: 

r 

The normalized eigenfunctions {vl^{r),u'^j^{r)} and the eigenvalues A^y (Table 1) are taken 
from our previous papers [22], [23]- Some lowest eigenfunctions {f Jy(r), m]y(?")} are plotted 



Mt = 0, r) = e yiAi^e-*-M/2 ^j^^^)^ ^^^^ = 0, r) = e ^/^— uiv(0 ■ (19) 



in Fig.3 (top), as they enter into the relevant combinations in the formula (18): 
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The background functions fN^r), $Ar(r) and their unstable eigenfunctions vlf{r),u^jy{r) have 
the following behavior: 



/^(r = 0) = l, /;(r = 0) = 0, /^(r^oo) = ±l, /^(r ^ oo) = 0, (21) 

$^(r = 0)=0, $'^(r = 0) = 0, <l>,v(r ^ oo) = <l;v, $'^(r ^ oo) = 0, (22) 

v'^ir = 0) = 0, v'^'{r = 0) = 0, v'j^{r ^ oo) = 0, v'j^\r ^ oo) = 0, (23) 

<(r = 0)=M*;v, u'n ir = 0) = 0, u'j^ir ^ oo) = 0, ^^'(r ^ oo) = 0, (24) 

where $Ar and the obtained constants. 

As a result, the boundary conditions for the considered evolutional Cauchy problem are 
summarized as follows: 

/(t,r = 0) = l, /,(t,r = 0) = 0, $(t = 0,r = 0) = SV, $,(t,r = 0) = 0; (25) 
lim/(t,r) = ±l, lim /,(t,r) = 0, 

r— >oo r— >oo 

Mm ^{t,r) = ^N, lim $,.(t,r) = 0. (26) 

r— >oo r— >oo 

Note that the value of the dilaton function at the origin is defined by the boundary 
conditions at the initial time only: $(t = 0,r = 0) = £t^. However, on each other slice 
t > this value is a free parameter that evolves according to the equations (|4]). In other 
words, we put free boundary condition for the dilaton function at the origin $(t > 0, r = 0). 
At the same time we keep the value of the dilaton fixed at infinity equal to its initial value 
limr_^oo ^{t > 0,r) = that provides the same scale during the evolution. 

Using the normalization condition for the unstable eigenfunctions {vlf{r),u^j^{r)} [22] 
one can obtain a useful formula for the mass-energy shift caused by adding eigenfunctions 



{f]v(r), M]y(r)} to the background static solutions in accordance with (18) 



E = ^^ + 2e2|Aiv|, (27) 

where is the mass-energy of the corresponding static solution {/jv(r), $jv(r)} and e is a 
free small parameter. As we will see below, the sign of the free small parameter e plays a key 
role for the late-time solution asymptotics. However, it enters squared into the mass-energy 



shift formula (27) in the framework of the considered perturbation scheme. 

So we have a well-posed Cauchy problem of mixed type (both, initial and boundary 
conditions are imposed), which has been solved numerically at various initial conditions 
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with the help of an adaptive mesh refinement of the finite difference scheme that preserves 
the total energy during the evolution (see P] for further details). The semi-infinite interval 
r G [0, oo) was replaced by an appropriate finite interval r G [0, -Rqo] for the numerical 
reasons. 



A. Decay of the A'^ = 1 static solution. 

Let us consider the = 1 static solution perturbed by its single unstable eigenmode 
with X\ = —9.0566 x 10^^ as initial data: 

/(t = 0,r) = /i(r) + ee-*^(")/2 i^^N = g, r) = $i(r) + e — ^^(r) , (28) 

r 



Mt = 0, r) = e yiAje-*^ v^^{r), $,(t = 0, r) = e «}(^) • (29) 

The results of various simulations with different values of e allow us to conclude that the 
late-time asymptotics of the solutions depends on the sign of the parameter e only. In fact, 
if the parameter e < is negative, the solutions evolve for a while in the linear regime 



(11), then start to collapse and later on they attain the stable N = self-similar attractor. 
Finally, the singularity develops at the origin in a finite time T. The absolute value of T 
depends on the value of e in the initial data. All such solutions with various initial e < 
have the universal self-similar late-time asymptotics with the effective dependence on the 
difference of the time t and the absolute blowup time T as r/{T — t). 

To recognize the self-similar nature of the dilaton function, one should extract its self- 



similar part (f){t,r) according to the (12), which corresponds to the frame appropriately 
shifted and rotated at the angle /3 = arctg{—2) back clockwise (see [9] for details). 

Typical profiles of the YM and the dilaton functions at various times are shown in Fig. 4. 
One should indicate the fact that the corresponding radial momentum density function ([s]) 
occurred to be always positively defined, starting with the initial time t = (see Fig. 3, 
bottom) until the blowup. This means that there is no any outgoing energy fiux at any time 
during the evolution. 



If e > the solutions evolve in the linear regime (11) for a while, then the nonlinear 
waves start moving towards a large r and finally one has a purely outgoing linear waves as a 
late-time solutions asymptotics. Indeed, as the outgoing wave propagates towards r — > +oo, 
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the influence of tlie nonlinear terms in the right hand sides of the equations (|4]) becomes 
negligible and one deals with the linear system of wave equations as follows: 



ftt + - frr - fr^r = 0, (30) 

2$ 

- ^ = 0. (31) 



The solutions of the linear system (30) have the standard form 



f{t,r) = fUt-r), <^(t,r) = ^^^^^^^, (32) 

where fiin{t — r) and (f)iin{t — r) are arbitrary functions. The asymptotical profiles of the 
relevant YM junit — r) and the dilaton (f)iin{t — r) functions are shown in Fig. 5 in the limit 
e — * 0. One should stress that similarly to the previous case e < 0, for e > the only 
purely outgoing energy flux was observed at any time of the evolution. In this respect the 
decay of the N = 1 static YMd solutions via their unstable eigenmode excitation is rather 
exceptional. 

B. Decay of the N = 2 static solution. 

Now we consider first the static N = 2 solution perturbed by its main unstable eigenmode 



with A2 = —7.5382 x 10 ^ and with the eigenfunctions {v\{r)^u\{r)} in the initial data (18), 



(119j). 

The obtained results occurred to be qualitatively the same as in the previously described 
case of the static = 1 solution decay. For e < the corresponding solutions start to 
collapse and later on they attain the = self-similar attractor that leads to finite-time 
blowup, see Fig. 6. However, in contrast to = 1 decay some outgoing energy flux is formed 
at the initial time as a sequence of the unstable eigenfunction shape at its outer side (see 
Fig. 3, bottom). One can observe that this outer part of the solution develops soon and then 
starts to propagate towards the spatial infinity separately on the main part of the solution 
(Fig. 6: bottom, right). 

For e > the solutions radiate away to the spatial infinity. Fig. 7 shows the typical 
pictures of such solutions. 

If we consider excitation of the second unstable eigenmode with \\ = —2.0742 x 10~^ and 



with the eigenfunctions {f|(r),U2(r)} in the initial data (18), (19), the late-time solutions 

14 



asymptotics occurred to be independent on the sign of the parameter e. Indeed, after some 
small oscillations near the static N — 2 configuration (see Fig . 8, 9) both solutions start to 
collapse, then they attain the N = self-similar attractor that leads to a finite-time blowup. 
In other words, the excitation of second unstable eigenmode in the N = 2 static solution 
does not provide a blowup/scattering alternative for the late-time solutions asymptotics. 
As we will show below, the same is true for all higher unstable eigenmodes of other static 
solutions with N > 1: the independent excitation of higher unstable eigenmodes always 
leads to the finite-time blowup. 

Note that small but detectable outgoing energy flux was observed in both cases e > 0, 
e < if the second unstable mode was excited. Similar to the main unstable mode excitation 
the outgoing pulse is formed already at the initial time owing to the unstable eigenfunctions 
shape at their outer sides. 

C. Decay of the AT = 3,4 static solutions. 

The static N — 3 and N — A solutions have 3 and 4 different unstable eigenmodes, 
respectively. 

An excitation of the main unstable modes with A3 = —4.9346 x 10^^ for the = 3 
YMd static solution and with A4 = —4.3455 x 10~^ for the A^ = 4 static solution provide a 
blowup/scattering alternative where the sign of e is a key parameter, similar to the A" = 1, 2 
static YMd solutions decay. 

Similarly to the static N — 2 solution decay, in the case of negative values of e, which 
leads to blowup, a small but detectable outgoing energy flux was observed arising at the 
initial moment of time, see Fig. 10. 

Along the same lines, the independent excitations of the higher unstable modes (with 
XI = -1.4957 X 10-^ Ai = -1.9622 x 10"^ for the AT = 3 static solution and with A| = 
-5.9905 x 10-^ XI = -1.3278 x 10"^ A^ = -1. x 10"^ for the A^ = 4 static solution) lead to 
the flnite-time blowup independently on the sign of the parameter e in the initial data. As 
a result, the asymptotical behavior of these solutions is independent on the sign of e, quite 
similar to the N = 2 static solutions decay, discussed previously. Similarly to the decay of 
the N = 2 static solution via their second unstable mode excitation, a small but detectable 
outgoing energy flux was observed in both cases e > 0, e < if the highest unstable modes 
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have been excited. 

We have also studied the possibihty of attaining some lower A^— 1, ..1 YMd static solutions 
following the decay of the static YMd solutions with zeros of the YM function (A^ = 2, 3, 4) 
in both ways (blowup/scattering). One can conclude that these solutions never attain other 
possible intermediate attractors of the low mass-energy. 

It seems reasonable that the decay of higher > 4 static YMd solutions is qualita- 
tively similar to the studied ones with A^ = 1,2,3,4. This allows us to draw the following 
conclusions. 



V. CONCLUSIONS AND DISCUSSION. 

We obtained strong numerical evidences in favour of the following statements: 

- All static YMd solutions are distinct local threshold configurations, separating blowup 
and scattering solutions. 

- The main unstable eigenmodes are only those responsible for the blowup/ scattering 
alternative in the final solution asymptotics for all static YMd solutions. 



- A sign of the small parameter e in (18) is a key parameter, providing blowup/scattering 
alternative when the main unstable eigenmodes are excited. 

- Excitation of higher unstable eigenmodes always leads to finite-time blowup indepen- 
dently on the sign of e. 

- Decay of the A^ = 1 static YMd solution via excitation of its unique unstable eigenmode 
is an exceptional case because the resulting waves revealing features peculiar to solitons. It 
was shown that such waves propagate in both radial directions as a whole without energy 
dispersion. 

Following the Bizon- Tabor conjecture [6j|, we believe that singularity formation key prop- 
erties are not specific features of definite nonlinear evolution equation, but rather are uni- 
versal characteristics of a criticality class which this equation belongs to. In this way the 
results obtained for blowup can be valid as well for all supercritical systems of the same 
type and can shed new light on the old problems inspired, in particular, by the last decade 
studies of black hole formation as a result of the collapse of massless matter fields. 

Let us remind that if gravity is included, say in the Einstein- Yang-Mills (EYM) sys- 
tem of spherically symmetric equations [2], a smallest black hole mass formed in Type-I 
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gravitational collapse occurs to be exactly equal to the mass of the lowest unstable = 1 
Bartnick-McKinnon (BK) solution. It is possible to take for granted that = 1 BK solution 
is a local threshold configuration which separates dispersing and collapsing solutions. Thus, 
the smallest black hole in Type-I collapse is formed if the unique unstable eigenmode of 
= 1 BK solution is excited with an appropriate sign of e in a limiting case as e — > 0. 
If these statements are accepted, then a natural question arises of why all matter which 
forms = 1 BK solution appears inside black hole event horizon and no any outgoing 
energy flux is produced. In fact, tales of outgoing pulses formed due to backscattering of 
ingoing collapsing waves on a space-time curvature are most general and universal sources 
of information about collapse details for a distant observer. In addition to this a various 
basic processes inside black holes, such as a mass inflation [21], [25] phenomenon, are also 
inspired by the backscattering of ingoing matter pulses. In the case under consideration an 
absence of such outgoing energy flux during A^ = 1 BK solution decay is very unusual and 
interesting. 

On the basis of the obtained conclusions, we believe that any supercritical system (with 
gravity included or not) has similar origin of this phenomenon. As we have shown, decay of 
the lowest static solution in YMd system via its single unstable eigenmode excitation pro- 
duces exceptional nonlinear soliton-like wave, which can propagate in both radial directions 
as a whole without backscattering of the energy. Accepting the universality conjecture it 
is possible to assume that similar soliton-like wave is formed in Einstein- Yang-Mills system 
that naturally explains the equality of the smallest black hole mass to the mass of the A^ = 1 
BK solution. We plan to study soliton-like behavior of such solutions in more detail in a 
separate paper. 

On the other hand, the collapse decay of higher A^ > 1 static unstable solutions is always 
accompanied by some energy flux to spatial infinity. Thus, the question about lower bound 
of mass-energy kink, attaining some distinct A^ > 1 static threshold configuration prior to 
blowup, is more subtle. We plan to consider this problem separately as well. 
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FIG. 1: Positions of YM functions f{t,r) (left) and the dilaton functions <I>(t, r) (right) with 
various initial data, shown in the moment of bouncing. Virtual positions of the = 1,2 static 
YMd solutions are also plotted in the same scale. One can see that the virtual = 1 static YMd 
solution plays a role of border line in space; all dispersing solutions can not penetrate this border 
line being bouncing outside. 
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FIG. 2: Evolution of the kink-type initial data for the YM function /(t, r) with the key parameter 
values a, closed to the threshold one: a < acnt. (left, top), a > acrit. (right, top). The mass-energy 
access which is approximately equal to the difference between the kink mass-energy and the mass- 
energy of the = 1 static solution, is radiated away towards r ^ oo before the solution attains this 
intermediate attractor. The plot of the corresponding radial momentum density function p{t, r) is 
also presented (bottom): the negative p{t,r) corresponds to the outgoing energy flux. 
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FIG. 3: The main unstable eigenfunctions entering in the relevant combination (20) for the static 



YMd solutions with = 1,2,3,4 (top). The radial momentum density function p{t = 0,r) at 



the initial moment t = in the Cauchy problem (18), (19) for the static YMd solutions with 
N = 1,2,3,4 perturbed by their main unstable eigenfunctions (bottom , left) with e < 0. The 
same is plotted (bottom, right) with the increased resolution in a vicinity of the right edge of the 
pulse. One can see that the only = 1 perturbed static YMd solution has a definite sign of radial 
density function p{t = 0, r) everywhere. 
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FIG. 4: The decay of the = 1 static solution along its single unstable eigenmode with Aj^ = 
—9.0566 X 10^^ and e = —0.2 which leads to blowup. The YM f{t, r) profiles (top, left) and the 
dilaton profiles ^{t,r) (top, right) at the various moments of time are plotted versus the logarithm 
of r. Virtual positions of the = 2, 3 static YMd solutions at the same scale are also shown. At 
the later time, prior to blowup, the solution attains the stable self-similar A^ = profile. The radial 
momentum density function p(t, r) (bottom, left) and the energy density function e(t, r) (bottom, 
right) are displayed at the same moments of time, which exhibits the purely ingoing energy flux. 
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FIG. 5: The asymptotical profiles of the YM function fiin{t, r) (left) and profiles of the normalized 
dilaton function (j)iin{t — r) = ($(t,r) — $oo) x f (^oo = 3.35879184 is the asymptotical value of 
the dilaton function for the N = 1 static solution) (right) are displayed at later moments of time 
versus retarded null coordinate r — t. These asymptotical profiles are the result of the N = 1 static 
solution decay along its single unstable eigenmode with A{ = —9.0566 x 10~^ in the limit e ^ 0. 
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FIG. 6: The decay of the N = 2 static solution along its main unstable eigenmode with A2 = 
-7.5382 X 10^2 and e = -0,2, which leads to blowup. The YM f{t,r) profiles (top, left) and 
the dilaton profiles <I>(t, r) (top, right) at the various moments of time are plotted versus the 
logarithm of r. At the later times, prior to blowup, the solution attains the stable self-similar 
= profile. The radial momentum density function p{t, r) (bottom: left, right) is displayed at 
the same moments of time. The negative values of the function p(t, r) (bottom, right) indicate an 
outgoing energy flux, which can be detected at spatial infinity. 



25 




FIG. 7: The decay of the N = 2 static solution along its main unstable eigenmode with A2 = 
-7.5382 X 10^2 and e = +0.1, which leads to scattering. The YM f{t,r) (top, left) and the 
dilaton $(t,r) (top, right) profiles are displayed at various moments of time versus the logarithm 
of the r. The radial momentum density function p{t, r) (middle, left) and the energy density 
function e(t, r) (middle, right) are shown at the same moments of time. The asymptotical profiles 
of the YM function fiinit,r) (bottom, left) and the profiles of the normalized dilaton function 
4'iin{'t — r) = {^{t, r) — $00) X T (^00 = 6.7478 is the asymptotical value of the dilaton function for 
the N = 2 static solution) (bottom, right) are also displayed versus the retarded null coordinate 
r-t. 
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FIG. 8: The decay of the N = 2 static solution along its second unstable eigenmode with = 
—7.5382 X 10^^ and e = ±0.5. The deviations from the static N = 2 solutions are shown for the YM 
f(t, r) — f]\f=2{r) (left) and for the dilaton <I>(i, r) — ^N=2{f') (right) functions at various moments 
of time at the beginning of evolution. Solid lines correspond to the solution with e = +0.5, dashed 
lines - to the solution with e = —0.5. 
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FIG. 9: The decay of the N = 2 static solution along its second unstable eigenmode with e = +0.5 
(solid lines) and e = —0.5 (dashed lines). The radial momentum density function p{t,r) (bottom, 
left) and the energy density function e{t, r) (bottom, right) are displayed at the same moments of 
time as in the previous Fig. 8. 
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FIG. 10: The decay of the A'^ = 3 static sokition along its main unstable eigenmode with A2 = 
-4.9346 X 10^2 and e = -0.2, which leads to blowup. The YM f{t,r) profiles (top, left) and 
the dilaton profiles ^{t, r) are plotted at various moments of time. The radial momentum density 
function p{t, r) (bottom: left, right) is displayed at the same moments of time. The negative values 
of function p{t, r) (bottom, right) indicate an outgoing energy flux, which can be detected at spatial 
infinity. 
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